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Abstract 

We present a comprehensive and self-contained discussion of the use of the transfer matrix to 
study propagation in one-dimensional lossless systems, including a variety of examples, such as 
superlattices, photonic crystals, and optical resonators. In all these cases, the transfer matrix has 
the same algebraic properties as the Lorentz group in a (2h- l)-dimensional spacetime, as well as 
the group of unimodular real matrices underlying the structure of the abed law, which explains 
many subtle details. We elaborate on the geometrical interpretation of the transfer-matrix action 
as a mapping on the unit disk and apply a simple trace criterion to classify the systems into 
three types with very different geometrical and physical properties. This approach is applied to 
some practical examples and, in particular, an alternative framework to deal with periodic (and 
quasiperiodic) systems is proposed. 

Keywords: Transfer matrix, hyperbolic geometry, periodic systems 



Contents 



1 Introduction H] 

2 Transfer matrix in quantum mechanics |3] 

2.1 Basic concepts on transfer matrix [3] 

2.2 Building the transfer matrix Q 

2.3 Hyperbolic Stokes parameters [8] 

2.4 Lorentz transformation associated to a transfer matrix [10] 

3 Transfer matrix in other contexts |12] 

3.1 Mechanical waves [12] 

3.2 Electromagnetic waves [12] 

3.3 Geometrical optics [14] 

4 The geometry of the transfer matrix [16] 

4. 1 Transfer function in the unit disc [16] 

4.2 Transfer function in the half -plane [20] 

4.3 Factoring the transfer matrix [22] 

4.4 Geometrical reflections as building blocks |23] 



Preprint submitted to Physics Reports 



today 



5 A closer look at the composition of transfer matrices |26] 

5.1 Setting up the inverse system |26] 

5.2 The Wigner angle |2i 

5.3 Hyperbolic turns |30] 

6 Periodic systems |32] 

6.1 Finite periodic structures |32] 

6.2 Bandgaps in the unit disc |34] 

6.3 Bandgaps in the half-plane |36] 

6.4 Quasiperiodic sequences |38] 

6.5 Hyperbolic tilings |42] 

7 Concluding remarks |44] 



1. Introduction 



Quantum mechanical scattering in one dimension describes many actual phenomena to a 
good approximation. The advantage of this theory is that it does not need special mathematical 
functions, while still retaining sufficient complexity to illustrate the pertinent physical concepts. 
It is therefore not surprising that there have been many articles dealing with various aspects of 
such scattering at various levels of detail, ranging from pedagogical issues ( |Eberly| 1965[ James 



1970! 'Formanek' TW; 'Kamal' TW; 'van Dijk and Kiers' 199T; Nogami and Ross', TW; 'Bar 
lette et al., 2000, 2001, Cattapan and Maglione, ,2003,, Sanchez-Soto et al. , 2005; Boonsem 



1983 



Jaworski and Wardlaw 1989 



Trzeci- 



and Visser[ |2010b) to edge-cutting research (Peres 

^owski and Gurioli 1993'; 'Sassoli-de-Bianchi, 1994'; 'Rozman et al. T994a'; 'Nockel and Stonej 
1994; Sassoli-de-Bianchi and Ventra 1995 ; Chebotarev and Tchebotareva , 1996 ; Kiers and v an] 



Dijk 



2000 



1996 Marinov and Segev 



Grosseletal. 2002 Boya 



1996^ 



Kerimov and Sezgin 1998 Visser 1999^ Miyazawa] 



2008t|Xuereb et al.||2009| |Boonser m and Visser||2009||2010; 



iwa] 
Oaf! 



These papers emphasize notions such as partial-wave decomposition, Lippmann-Schwinger 
integral equations, transition operator, or parity-eigenstate representation, paralleling as much as 
possible their analogues in two and three dimensions. In other words, these approaches, like most 
of the standard textboo ks on the subject ( . Goldberg er and Watson 1964, Newton^ 1966, CoherTl 



Tannoudji et al. 1977 Galindo and Pascual 1990[l, employ the S matrix [note some significant 



exceptions, such as, e.g. Mathews and Venkatesan 
or[S5ghl([T997l)]. 



1978 1, Merzbacher ( 1997 1, Ballentine ( 1998 1, 



The elegance and power of the S-matrix formulation is beyond doubt. However, it is a 
"black-box" theory: the system under study is isolated and is tested through asymptotic states. 
This is well suited for experiments in elementary particle physics, but becomes inadequate as 
soon as one couples several systems. The most effective technique for studying such a coupling 
is the transfer matrix, in which the amplitudes of two fundamental solutions on either sides of a 
potential are connected by a matrix M. 

The transfer matrix is a fruitful object widely used in the treatment of layered systems, like 

fT98F, 'Ram-Mohan et al. 



superlattices (Tsu and Esaki 1973 Esaki 



iiengl [T9891 |Vinter and Weisbuch 1991 Weber 1994; Sprung et al. 



19881 Hauge and Stov- 



2003 I or photonic crys- 



tals piaannopoulos et al. 1995 Bendickson et al. 1996 |Tsai et al. 1998 1. Optics, of course, is 



a field in which multilayers are important and the method is time honored ( Brekovskikh 
|Lekner|[T987i|Azzam and Bashara|[T987l|^[T988l l. 
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1960 



An extensive and up-to-date review of the applications of the transfer matrix to many prob- 



lems can be found in the two excellent m onographs by Garcia-Moliner and Velasco ( 1992 1 and 



Perez-Alvarez and Garcia-Moliner (2004 1. They are addressed to anyone who wants to enter the 
field and provide a really professional level of penetration into the basic issues. 

A natural question thus arises: why yet another essay on the transfer matrix? The answer is 
simple: a quick look at the literature immediately reveals the different backgrounds and habits 
in which the transfer matrix is used and the very little "cross talk" between them. In fact, many 
scientists are usually not aware of the mathematical basis behind the standard toolkits they are 
using in their everyday research. The main goal of this review is precisely to fill this gap. 

When one thinks in a unifying mathematical scenario, geometry immediately comes to mind. 
Although special relativity is the archetypal example of the interplay between physics and ge- 
ometry, one cannot forget that geometrical ideas are essential in the development of modern 
physics ( |Schutz| [1997] [Kaudererj [1994^. 

In recent years a number of geometrical concepts have been exploited to gain further insights 



into the behavior of scattering in one dimension ( 


Yonte et al. 


2002 


Monzon et al. 


2002 Barriuso 


etal. 2003a 


2004 ; Sprung et al. 2004 Martorel 


I etal. 2004 


Sane 


lez-Soto et al. 


2005 ; Barriuso 


et al. 2009 


1. The algebraic basis for these developments is the fact that the transfer matrix 



is an element of the group SU(1, 1), which is locally isomorphic to the {2 + l)-dimensional 
Lorentz group SO(2, 1). This leads to a natural and complete identification between reflection 
and transmission coefficients and the parameters of the corresponding Lorentz transformation. 



As soon as one realizes that SU(1, 1) is also the basic group of hyperbolic geometry (Coxeter 
|1968'), it is tempting to look for an enriching geometrical interpretation. In fact, we propose to 
look at the the action of the transfer matrix as a bilinear (or Mobius) transformation on the unit 
disk, obtained by stereographic projection of the unit hyperboloid associated with SO(2, 1). 

Borrowing elementary techniques of hyperbolic geometry, we can classify and reinterpret all 
the relevant features of these matrices in a very elegant and concise way, largely independent 
of the model considered. We stress that this formulation does not offer any inherent advantage 
in terms of efficiency in solving practical problems; rather, we expect that it could supply a 
general and unifying setting to analyze the transfer matrix in many fields of physics, which, in 
our opinion, is more than a curiosity. 



2. Transfer matrix in quantum meclianics 



2.1. Basic concepts on transfer matrix 

We consider the quantum scattering of a particle of mass m in one spatial dimension by a 
potential barrier V{x). This is governed by the time-independent Schrodinger equation 



where 



2m _ 



>F(x) = e1'(x), 



U(X): 



2m 



-V(x), 



(1) 



(2) 



E being the energy of the particle. We assume this potential to be real (i. e., a Hermitian operator) 
but otherwise arbitrary in a finite interval (a,b). Outside this interval it is taken to be a constant 
that we define to be the zero of energy. Complex potentials can be used to model absorption, a 
situation which we shall not touch upon ( |Monz6n et al.)|201 
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The treatment can be also adapted, with minor modifications, to deal with potentials for 
which Va + Vh, and also with the more subtle case of a pos ition-depe ndent effective mass m(x), 
which usually arise in superlattices (Leibler,,1975 Bastard 1981;, Perez- Alvarez and Rodriguez 



Coppola! [T988||Thomsen et al.[[T989[]Burt||1992| i. 

Since £ > 0, the spectrum is continuous and we have two linearly independent solutions for 
given value of E ( Galindo and Pascuall|1990| . Accordingly, the general solution of the time 



independent Schrodinger equation can be expressed as a superposition of a right-mover e 
a left-mover e~'*^^: 



+ik.\ 



^ +ik(x-a) ^j^_g-ik{x-a) 



^ab(x) a < X < b 

g +ik(x-h) ^ B_e-ik(x-b) ^ > ^ 



and 



(3) 



where - s and the subscripts -i- and - indicate that the waves propagate to the right and to 
the left, respectively (see figure 1). The origins of the movers have been chosen to simplify the 
subsequent calculations. 

To solve the problem in a closed form one must work out the Schrodinger equation in (a,b) 
to compute ^abix) and invoke the appropriate boundary conditions, involving not only the conti- 
nuity of itself, but also of its derivative. In this way, one has two linear relations among the 
coefficients A+ and B+, which can be solved for any amplitude pair in terms of the other two: the 
result can be expressed as a matrix equation, which translates the linearity of the problem. For 
our purposes, it is more advantageous to express a linear relation between the wave amplitudes 
on both sides of the scatterer, namely. 



ah 



B- 



(4) 



Maf, being the transfer matrix for the potential. 

The reflection and transmission coefficients are the ratio of the amplitudes of the reflected and 
transmitted waves to the amplitude of the incoming wave, respectively. Denoting the amplitudes 
for waves propagating from the right as r/^, and tha and repeating the procedure, one easily finds 
that time-reversal invariance imposes 



^ab — ^ba-! ^ba/^ba — ^abl^ab^ 
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(5) 



while the conservation of the flux (??) gives 



k„/,p + luP = i. (6) 

In conclusion, the form of the transfer matrix is 

In the particular case of a symmetric potential [i.e., V{x) = V(-x)], it is clear that rah = rba and 
therefore the matrix element /3 is an imaginary number. 

Thus far, we have related the amplitudes A± to the B±, as in equation This choice is by no 
means essential and we could relate the amplitudes taken in the reverse order. The corresponding 
transfer matrix, represented by M^a, can be expressed as 

Mba-l ?t (8) 



ab 



where we have used (|5]). 

We will now bring up the paradigmatic example when the potential V{x) reduces to a rectan- 
gular potential barrier of width L and height Vq. The calculations can be easily carried out, so we 
skip the details (Cohen-Tannoudji et al. 1977[ l and simply quote the results for rab and t^i, (with 
the choice of movers in figure [TJ 

rab 



(k?- + K^) sin(A:L) + 2ikKC0s{KL) ' 

(9) 

2ikK 

(k^ + K^) sm(KL) + 2ikKCOsh{KL) ' 

with - 2m(E-Vo)/h^- These coefficients correspond to E >Vq. When E <Vq the expressions 
are 

(A;2+/f2)sinh(^L) 



rab = 



tab - 



(k^ - K^) sinh(/cL) + 2i^/ccosh(/fL) ' 

(10) 

2ikK 



(k^ - K^) sinh(/cL) + 2/^/ccosh(/cL) ' 



where now k - 2m(VQ-E)/h^. This can be obtained from the previous case with the formal 
substitution k — > k. Finally, when E = Vo a limiting procedure yields 

kL 2i 
''"'~kL + 2i' '"''~kL + 2i- ^ ' 



A detailed discussion of the significance of these three situations can be found in Bohm ( 1989 1. 

It is worth stressing that one could also relate outgoing amplitudes in terms of the incoming 
amplitudes (which are the magnitudes one can externally control). This is precisely the scattering 
matrix, which can be concisely written as 

^ Sabi ^ ], (12) 



A_ / """l B 
5 



and Sab reads 



^ ab ^ba 



(13) 



Due to the properties (|5]l and (j6]), Sab is unitary and with unit determinant, that is, an element of 
the group SU(2). 



N ote carefully that the transfer matrix depends on the choice of basis vectors ( Perez- Alvarez 
et al. 2001 Perez- Alvarez and Garcia-Moliner 2004 \ and special care must be paid when com- 
paring results from different sources. For example, instead of specifying the amplitudes of the 
right and left-moving waves, we could also write a linear relation between the values of the 



wave function and its derivative at the points a and b ( Sprung et al. 1993 i. prefer to employ the 
adimensional variables 



*P(x) + -r(jc), 
k 



k 



(14) 



This looks very much as passing from position-momentum to creation-annihilation operators ( Bal 



lentine 1998 1. The amplitudes associated to these variables are related to A± by 



with 



04- 



1 



and analogously at the point b. The transfer matrix in this equivalent representation is 

Mab^'UUab'U^ 



(15) 



(16) 



(17) 



where t stands for the Hermitian conjugate and 

a = ReaH-Im/?, 



(18) 



c = -Imc + Re/?, h -R&a -Imfi , 



are real numbers. The role that the transformation U will play in what follows justifies our 
choice in ( [T4] i. 

Since the determinant are preserved by matrix conjugation, we have that detAlu/, - + \. In 
other words, the matrices Mab belong to the group SL(2, R) of unimodular 2x2 matrices with 
real elements. The transformation by U establishes in fact a one-to-one map between the group 
SL(2, R) of matrices Mat and the group SU(1, 1) of matrices Mab, which allows for a direct 
translation of the properties from one to the other, as we will have occasion to check. 

Irrespective of the representation used, transfer matrices are very convenient mathematical 
objects. Suppose we know how the wave functions "propagate" from point b to point a, with a 
transfer matrix we symbolically write as Mn^, and also from c to b, with Mfof The crucial point 
is that the propagation from c to a is described by the product 



Uac = Mfli Ubc . 
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(19) 



This property is rather helpful: we can connect simple scatterers to create an intricate potential 



landscape and determine its transfer matrix by simple multiplication (Jonsson and Eng [1990 



Kalotas andTee] [T99T] [Walker and Gathrightj [T994} [Rozman et aL] |1994a[|Yuan et al.| |2010f . 



However, this important property does not seem to carry over into the scattering matrix in any 



simple way (Aktosun 1992 Aktosun et al. 1996 Giust et al. 2009 because the incoming 
amplitudes for the overall system cannot be obtained in terms of the incoming amplitudes for 
every subsystem. While this is not a difficulty for a single scatterer (the typical situation arising 
in particle physics and for which the S-approach was specially tailored), it constitutes a drawback 
in applications where a number of cascaded systems are present. 



2.2. Building the transfer matrix 

The complete determination of the transfer matrix for an arbitrary potential Y(x) amounts 
to solving the Schrodinger equation and, in consequence, it is not, in general, a simple exercise. 



1995 



et al 



Very accurate approximation schemes are available, among which we cite the WKB dChebotarev 
1997), the variational ( [Bastard et al] [T983| |Ahn and Chuang] [1986} pouiai ]l 995 HAndo 



2003 I, the Monte-Carlo (Singh 1986 ; Kalos and Whitlock"2007), or the finite-element 



methods (Hayata et al. 1988 Nakamma et al.,, 1989,, Ram-Mohan,, 2002, |Liu et al. 2004 1. 



Here, we favor a method inspired by the long experience in dealing with layered systems ( Ken 


nett 1983 ; Perez-Alvarez et al. 1988 ; Rodriguez-Coppola et al. 


1990 Garcia-Moliner and Ve- 


lasco 1992 


,Perez-Alvarez and Garcia-Moliner, ,2004,). Roughly speaking, the idea is that one 



can consider V(x) as made of successive constant barriers, as schematized in figure |2] ( Kalotas[ 



and Lee 1991 ; Rozman et al. 



He et al.l 



2 005,.Monsoriu et a 



1994b; Grossel et al. , 1994; Cao et al. 2001 ; Rakityansky , 2004 
.l ,2005,.Su et al.,,2008.,Hutem and Sricheewin, .2008,. Wen et a~ 



20101. The yth barrier, of height V/ and width dj, is situated between the points Xj-\ and Xj (we 
take X{s-a and jca^+i - h). In this way, we express \Aab as a product of matrices that characterize 
the effects of the individual discontinuities and propagations of the entire discretized structure, 
taken in the proper order, as follows ( Yeh 1988[ l: 



ab = loiPlll2P2l23---lo-l)7PjliO+ 



+ 1) ■ ■ ■ I(A'-1)A'Pa'U(A'+1) ■ 



(20) 



When the number of barriers is large enough, the method should provide satisfactory results for 
any potential (Jirauschek 2009|l. Here l,j account s for the discontinuity at the interface between 
y, and Y j, and has the form (Landau and Lifshitz 200 l| l 



1 / 1 



1 



(21) 



where r,y and r,y are the reflection and transmission coefficients at the interface ij and are given 
by 

K; - Kj 2Ki 

nj - 



Ki +K 



K,+Ki 



(22) 



The wave number is Kj - 2m(E - Vj)/h and, for simplicity, we have assumed E > Vj. They 
verify 



detl,j - 



(23) 



and the outstanding composition law 



(24) 



V. 

-J 



^ 



V(x) 

A 

♦ 



'Kv^. 



a 



Figure 2: Decomposition of the potential V{x) in elementary rectangular barriers. 



The matrix Pj describes the effect of propagation through the barrier j alone, and reads 



sxpi-i6j) 
where the phase shift is dj - Kjdj. Now, we have 

det = 1 . (26) 



By taking the determinant in equation ( 20 1 and using ( 23 1 and ( [26] l we get a simple but 
relevant result: 

detM„/_, = = — . (27) 

KQ tab 

Therefore, when Va and Vb are the same (as we have assumed until now), the determinant of VAab 
is +1. When these potentials are different, this result also holds by renormahzing conveniently 
the amplitudes ("Monzon and Sanchez-Soto 1999[ l. 
If we denote 

M; = loiPil;0, (28) 

which corresponds to the jth barrier sandwiched between two identical constant potentials (that 
for simplicity we take as 0), we can also rewrite equation ( 20 1 as 

N 

Uab^W^j. (29) 

;=i 

Sometimes, it is more convenient to work with unimodular matrices in SL(2, R). The coun- 
terparts of interface and propagation can be obtained by conjugating pT| ) and ( [25] l with 1/; the 
final result is 

1/ 1 '"'./■ \ 03 . -i COS 5/ sin5j 



^ tii \ nj I I M -sin*^; cos 5; / 



2.3. Hyperbolic Stokes parameters 

To move ahead let us construct the matrices 



where X = A or B aie the amplitudes that determine the behavior at the points a and b, respec- 
tively. They are quite reminiscent of the coherence matrix in optics or the density matrix in 
quantum mechanics (Mandel and Wolf 1995| l. Observe that J is Hermitian and detJ - 0. In 
addition, one can readily verify that 

Ja^Mat^bMl^, (32) 

so they transform under Mat by congruence. 

Let now cr^ (the Greek indices run from to 3) be the set of four Hermitian matrices ctq-IL 
(the identity) and {0-1,0-2,0-^) (the standard Pauli matrices). They constitute a natural basis of 
the vector space of 2 x 2 complex matrices, so the coordinates 5'^ with respect to that basis are 



/ = -Tr(Ja-^), (33) 



so that 





= ^(|X+|2 + |X_|2), 

= Re(X;x_), 

= im(x;x_), 

= ^(|X+|2-|X_|2). (34) 



The congruence ( 32 1 induces in this manner a transformation on the variables of the form 

^=At^;, (35) 

where A^y can be found to be 

A''yiMab) = \ Tr{cr^Mab<TyMli^) , (36) 
and it turns out to be a Lorentz transformation. This equation can be solved to obtain Mai, from A. 



The matrices Mat and -M„i, generate the same A, so this homomorphism is two-to-one (Barut 
[T980| . 

The variables are coordinates in a Minkovskian space. Since det J = 0, the value of the 
interval (in both points a and b) is 

(,0)2 _ (^1)2 _ (^2)2 _ (^3^2 (37) 

SO it is lightlike. Moreover, the conservation of the probability current expressed in equation (??) 
means that the coordinate s^, defined in ( 34 1, remains invariant and ( 37 1 reduces to 

(i'')2-(i')2-(5Y = (i^)^ = constant, (38) 

that is, a two-sheeted hyperboloid of radius s^, which without loss of generality will be taken 
henceforth as unity (see figure [3]l. All this shows that the group SU(1,1) of transfer matrices is 
locally isomorphic to the (2+l)-dimensional Lorentz group SO(2,l). In fact, in more technical 
terms, the matrices Mab form a two-dimensional spinor representation of the restricted Lorentz 
group SO(2,l): A1A2 < — > M(Ai)M(A2) = +M(AiA2). 




Figure 3: Pseudospherical coordinates of tlie unit two-sheeted hyperboloid associated witli a transfer matrix M^, 



Let us now rewrite the complex amplitudes Xj. in polar form 

X± = |X+|expOV±). 



(39) 



Denoting (f - (p+ - (p- and introducing the angle x in pseudospherical coordinates (shown in 
figure [3]), we have 



= cosh;^' , 

- sinhx cos (p , 

- sinhx sin (p . 



(40) 



This parametrization is very evocative of the standard one for the Stokes parameters in the unit 
Poincare sphere ^Born and Wolf 1999 1, except for the fact that now the angle appears as the 
argument of hyperbolic functions. We refer to these parameters as hyperbolic Stokes parameters. 
One can interpret the rotations (p and the hyperbolic rotations x much in the same way as it 
is done for the Poincare sphere ( |Giust and Vigoureux 2002| l: a rotation of angle tp about the 
axis corresponds to a dephasing between left- and right-traveling waves without changing 
their relative amplitudes; as it happens, for example, in the propagation inside a barrier. On the 
contrary, a hyperbolic rotation of angle x corresponds to a change in amplitude between left- and 
right-traveling waves, as it occurs at the discontinuity between two barriers. 



2.4. Lorentz transformation associated to a transfer matrix 

To better appreciate the physical meaning of the Lorentz transformation induced by a transfer 
matrix, we recall that a Lorentz transformation A can be always decomposed into a product of a 
spatial rotation R and a boost L along an arbitrary direction ( |Moretti, ,2006j 

A = LR. (41) 

The analogous factorization for a transfer matrix is the polar decomposition, which ensures that 
any matrix M e SU(1,1) (to simplify the notation, we drop the subscript ab in the rest of this 
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section) can be expressed in a unique way as 

M = HU, 



(42) 



where H is positive definite Hermitian ("modulus") and U is unitary ("argument"). Under the 
homomorphism discussed previously, H generates a boost and U a rotation, in agreement with 



equation (41 



To find their explicit form we note that polar decomposition for the matrix M in equation (|7]l 
reads 



M 



HU 



1 



exp(-/T) 
exp(/T) 



where we have expressed the reflection and transmission coefficients in the form 

r - \r\ exp(/p) , t = \t\ exp(/T) . 

Using p6l ), the unitary component U generates the rotation [in (2+1) dimensions] 



R(U) = 



( 1 
cos(2t) 
sin(2T) 



^ 
- sin(2T) 
cos(2t) 



(43) 



(44) 



(45) 



that is, a spatial rotation in the plane 1-2 of angle twice the phase of the transmission coefficient. 
The Hermitian component H generates the boost 



L(H): 



y —yucosp —yvsmp 

-yvcosp l+(y-l)cos^p (y- l)cospsinp 
-yysinp (y- l)cospsinp l+(y-l)sin^p , 



(46) 



The modulus of the velocity v (we take c-1 everywhere) and the relativistic factor y = 1 / 
of this boost are 

2\r\ l+\rf 



V = 



l+|r| 



2 ' 



l-|r| 



2 ■ 



(47) 



The matrix L(H) is then a boost to a reference frame moving with a constant velocity u in the 
plane 1 - 2, in a direction forming a counterclockwise angle p with the axis 1. 
If, as it is usual, we introduce the rapidity ^ from ( Jackson 1975\ 



V = tanh^, 



(48) 



we have the following appealing identification of the reflection and transmission coefficients with 
the parameters of the Lorentz transformation: 



r = tanh(^/2)exp(/p). 



f - sech(^/2)exp(/T). 



(49) 



Therefore, \r\ - tanh(^/2), behaves as a velocity, while \t\ behaves as 1/y. 

The convenient properties of the hyperbolic tangent have been exploited in dealing with 
layered systems ( |Khashan]|1979[[Corzine et aL||1991j l, and explain why the reflection coefficients 
are examined in greater detail than the transmission ones. 
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3. Ti^ansfer matrix in other contexts 



One-dimensional continuous models provide a detailed account of the behavior of a variety 
of systems (|Lieb and Matthis 1966 Albeverio et al. 2004| l. The nature of the actual particles, or 
states, or elementary excitations, as they may be variously called, is irrelevant for many purposes: 
there is always two input and two output channels related by a 2 x 2 transfer matrix. In fact, this 
matrix can be viewed as a compact way of setting out the integration of the differential equations 
involved in the model with the pertinent boundary conditions; this is what makes the method so 
effective. 

From this perspective, one can construct a general theory of the transfer matrix for second- 
order differential equations (Khorasani and Adibi (2003 ; Khorasani and Mehrany 2003 ). How 



ever, we prefer to explore some selected examples (Perez-Alvarez and Garcia-Moliner 2004 1; 
this restricted choice reflects the authors personal bias, but it is illustrative enough to grasp what 
the method is about. 



3.1. Mechanical waves 

Transverse waves on weighted strings, longitudinal waves on loaded rods, acoustic waves 



in corrugated tubes, and water waves crossing sandbars, among other examples (Griffiths and 



Steinke 2001 1, are governed by the classical wave equation 



(50) 



Here il/(x, t) is the amplitude of the considered phenomena (in the above-mentioned examples 
^{x, t) stands for the transverse displacement of the string, the displacement of a point whose 
equilibrium position is x, the pressure above ambient, or the height of the surface above its 
equilibrium level, respectively), and v(x) is the local propagation speed of the perturbation. We 
continue to use complex notation, with the understanding that the physical wave is the real part. 
For a monochromatic perturbation of angular frequency co [i/'(x,f) = 'V{x)e~'^^~\ equation (50 1 
reduces to 

-+k^{x) 



dx^ 



and the local wave number is 



k{x)^ 



'¥(x)^Q, 



v{x) 



(51) 



(52) 



The expansion in left- and right-movers in section 2 can be transplanted here without modifi- 
cations. In addition, as we did in 2.2, we can build the transfer matrix by assuming that the 
material parameters of the medium are piecewise constant and vary in a stepwise manner, so 
each constituent slab is a homogeneous material by itself. The mismatched impedances gen- 
erate the reflected and transmitted waves ( [Crawford} |1968| l, while the application of the proper 
boundary conditions at the discontinuity points (which depend on the particular model under 
consideration) provide the corresponding amplitude coefficients. 



3.2. Electromagnetic waves 

We next consider the propagation of plane electromagnetic waves in a (nonmagnetic) strati- 
fied medium, whose optical properties are contained in the dielectric function e{x) - n^(x) [n(x) 
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Substrate 



Figure 4: Amplitudes of the input [A+ and B_] and output [A_ and B+] fields in a multilayer sandwiched between two 
semi-infinite ambient and substrate identical media. The angle of refraction in the jth medium is denoted 0j. 



is the local value of the refractive index] . For a monochromatic component of frequency w we 



write the field components as (Monsivais et al. 1995 i 

E(r, f ) = &ix) exp[-/(wf - K ■ r)] , B(r, f) = S(x) exp[-/(wf - K ■ r)] , 



(53) 



where K is the component of the wave vector in the plane perpendicular to the x axis. By 
eliminating, e.g., the magnetic field B from Maxwell equations, it turns out that ( |Bom and Wolfj 
[19991 ) 

— +k^(x) 



dx^ 



£(x) = 0, 



and k{x) is the local value of the normal component of the wave vector 



k(x)^ 



e{x)%^-K\ 



(54) 



(55) 



A completely analogous equation can be written for the magnetic field by eliminating E. 



Apparently, equation ( 54 1 is identical to ( 5 1 1, and the theory can be immediately extended 
here, expressing the solution as a superposition of a left- and right-mover fields. But this requires 
some extra care because the amplitude in ( [54| is a vector 

For linear isotropic media, any plane wave can be written as a superposition of an s (or 
TE) wave and a p (or TM) wave. The s wave has its electric vector perpendicular to the plane of 
incidence, and the p wave has its electric vector in the plane of incidence (and its magnetic vector 
perpendicular to the plane of incidence; hence its designation as a TM, or transverse magnetic. 



wave). If we further take the plane of incidence to be the (x,z) plane, the vectors are ( Azzam and 
IBasharal [T987i |Lekner[ [T9871 [T988] l 








{s polarization) 



^ 





{p polarization) . 



(56) 



This has to be taken into account when matching the boundary conditions between two media. 
For these two basic polarizations the problem reduces to a scalar one, and the transfer matrix can 
be applied as before. 

If the medium extends from x-a\ax-h, bounded by the same homogeneous media (ambi- 
ent and substrate) of index «o, the formal solution can be written in full analogy with equation (|3]l 
and one can use a transfer matrix that relates the field amphtudes A± and B±. 

13 



Again, one can take the medium as consisting of a stack of 1,2,..., j,...,m plane-parallel 
layers, as sketched in figure|4] We denote by rij, dj, and 9j, respectively, the refractive index, the 
thickness, and the angle of refraction of the jth medium, which can be obtained by a repeated 
application of Snell's law (which is a consequence of the conservation of the modulus of K) 



no sin 00 = ■ ■ ■ nj sin 6j = ■■■ = «,„ sin 0„, . (57) 
The transfer matrix is given by the ordered product in equation ([20| and the interface matrix 



has the same expression as in equation (21 1, but with 



p rij cos Hi - rii cos tij 
'j n j cos Oj + Hi cos 0j ' 



n, cos Oi - rij cos 9j 
rtjcosdj + njCOS0j' 



Irij cos Oj 



'j rijCOsOi + rtiCOsOj' 
Irii cos 6i 



(58) 



'J 



rii cos 6i + rij cos 9j ' 



for each one of the basic polarizations. The propagation matrix is also as in (25 1, but now the 
phase shift is 

In 



— rtjdjCosGj, 



(59) 



A being the wavelength in vacuum. All this gives the formalism developed by Hayfield and White 
in terms of movers ( Azzam and Basharal|1987| l. 

It is also possible to develop an equivalent formalism by employing the amplitudes in equa- 
tion ( fTS) , which, roughly speaking, are the electric field and its derivative at each point. This is 
the idea behind the pioneering work of Abeles ( 1948[l. 



3.3. Geometrical optics 

Finally, we look at the paraxial propagation of light through axially symmetric systems, con- 
taining no tilted or misaligned elements ( |Wolf 2004 1. We take a Cartesian coordinate system 
whose X axis is along the axis of the optical system (see figure |5]l and represent a ray at a plane 
X by the transverse position vector q{x) (which can be chosen in the meridional plane) and by 
the momentum p(x) = n9(x), which in the paraxial limit is p(x) - ndq(x)/dx, where n is the re- 
fractive index of the medium. These are canonical coordinates and satisfy all the mathematical 
requirements for a consistent description (Guillemin and Sternberg 1984| ). 




I ^ > X 

Figure 5: Notation for the ray vector. The optical axis is x, q{x) is the transverse position at a point x reached by a ray 
and 6 is the ray inclination at that point. 
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Figure 6: Optical cavity consisting of two spherical mirrors of radii Ri and i?2 separated a distance d. 



In homogeneous media, to fully specify the ray behavior three basic matrices are needed, 
namely 

I d/n\ I I \ / 1 

Oil' \ in'-n)/R 1 I' 1 -2n/R 1 



(60) 



The first one gives the propagation through a distance d, the second gives the changes in the 
ray parameters for a refraction in a dioptre of radius R separating two homogeneous media of 
refractive indices n and n', and the third one is the reflection in a mirror of radius R. 

Let us apply these matrices to the illustrative example of an optical cavity consisting of two 
spherical mirrors of radii Ri and R2, separated a distance d, which will be examined in more 
detail in Section 6.3 (see figure 6| . The ray-transfe r matrix corresponding to a round trip can be 
routinely computed using (|60]) (Gerrard and Burch]|1975[): 



d \ 
2gig2-gi+g2-i ::(2gig2+gi+g2) 

2 

, -j('2.glg2-gl-g2) 2gig2+gl-g2-l 

V d 



(61) 



where we have introduced the parameters (i - 1 , 2) 

gi^l-- 



(62) 



In the same vein, a general first-order system can be built as a cascaded application of these 
three basic elements and the ray parameters change according to the simple transformation ( Si- 
mon and Wolf| [20001 [Woif| |2054t [Baikal et al.[ [2004) 1 



Pa \ Pb 



and Mab is the ray-transfer matrix 



Mab ■■ 



a b 
c b 



(63) 



(64) 



Since the three basic matrices in ( [60| l have unit determinant, this means that det A1 = +1, so that 
they belong to the group SL(2, R) of real unimodular 2x2 matrices. This is the essence of the 
celebrated abed law in geometrical optics. 
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4. The geometry of the transfer matrix 



4.1. Transfer function in the unit disc 

Let us go back to the unit two-sheeted hyperboloid ( [38| that is the phase space for our prob- 
lem. If one uses stereographic projection taking the south pole S = (-1,0,0) as projection center 
(see figure|7]i, the projection of the point (s", s^ , s^) becomes in the complex plane 

foT X - A or B. This confirms that what matters here are the transformation properties of ampli- 
tude quotients rather than the amplitudes themselves. In terms of the pseudospherical coordinates 



(40 1, this point can be written as 



sinh y 

Z = — exp(iV) = tanhOif /2) exp(/93) , (66) 

1 + cosh;^' 

which allows to interpret as a rapidity and (/? as a phase shift. The upper sheet of the unit 
hyperboloid is projected into the unit disc, we shall denote D, the lower sheet into the external 
region, while the infinity goes to the boundary of the unit disc. We mention in passing that 
stereographic projection is conformal, meaning that it preserves the angles at which curves cross 
each other on the two-sheeted hyperboloid. 

Through stereographic projection, the standard Minkowski distance in the unit hyperboloid 
becomes in D ( |Andersonl 1 1 999[ l 

, 7 dzdz* 

ds^ = . (67) 

(i-kP)2 

The geodesies in the hyperboloid are intersections with the hyperboloid of planes passing through 
the origin. Consequently, hyperbolic lines are obtained from these by stereographic projection 
and they correspond to circle arcs that orthogonally cut the boundary of the unit disk (diame- 



ters are a particular instance of these geodesies), as equation ( |67] i confirms after some calcula- 
tions (Misc henko and Fomenko|[l98 8). 



It seems natural to consider the complex variables in equation ( 65 i for both points a and b. 
The basic linear relation expressed in equation Q settles a transformation on the complex plane 
C, mapping the point Zh into the point Za according to 

Za = 0[M„fo,zfc] = , (68) 

pzb + a 

which is a bilinear or Mobius transformation. The action of the transfer matrix appears then as 
a function Zu - f{zb) that can be appropriately called the transfer function ( Yonte et al. 2002[ l. 



One can check that the unit disk D, the external region and the boundary remain also invariant 
under this action. 

We shall need the concept of hyperbolic distance in the unit disc. To this end, it is customary 
to define the cross ratio of four distinct points za, zb, zc, and zd as the number 

, I s (za-zc)/(zb-zc) 

(za,zb\zc,zd) = 7-. r , (69) 

{za-zd)/(zb-zd) 

which is real only when the four points lie on a circle or a straight line. In fact, bilinear transfor- 
mations preserve this cross ratio ( (Pedoe 1970| l. 
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Figure 7: Outline of the unit hyperboloid and a geodesic on it. We also show how a hyperbolic line is obtained in the 
unit disk by stereographic projection, taking the south pole as projection center. 



Let now z and z' be two points that are joined by the hyperbolic line whose endpoints on the 
unit circle are e and e' . The hyperbohc distance between z and z' is 



dii{z,z')^]-\\n{e,e'\z,z')\ 



(70) 



The fundamental point for us is that bilinear transformations are isometrics; i.e., they preserve 
this distance. 

The visual import of the disk with this metric is that a pair of points with a given distance 
between will appear to be closer and closer as their location approaches the boundary circle. Or, 
equivalently, a pair of points near the unit circle are actually farther apart (via the metric) that a 
pear near the center of the disc, which appear to be the same distance apart. 

An important tool for the classification of the transfer-matrix action are the fixed points, 
which correspond to the wave functions such that Za - Zh in equation ( [68] l. If we denote them by 
z/ we have that 

z/ = <D[M„i,z/], (71) 



whose solutions are 



1 



Zf± = — |-2/lm(a) ± ^[Tr(M„,)]2-4 



(72) 



When |Tr(M„i,)| < 2 the action is said elliptic and it has only one fixed point inside D, while 
the other lies outside. Since in the Euclidean geometry a rotation is characterized for having only 
one invariant point, this action can be appropriately called a hyperbolic rotation. 

When I Tr(M„f,)| > 2 the action is hyperbolic and it has two fixed points both on the boundary 
of D. The hyperbolic line joining these two fixed points remains invariant and thus, by analogy 
with the Euclidean case, this action will be named a hyperbolic translation. 

Finally, when | Tr(Maf,)| - 2 the action is parabolic and it has only one (double) fixed point on 
the boundary of D. This action has no Euclidean analogy and will be called a parallel displace- 
ment for reasons that will become clear soon. 
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In is worth mentioning that for the example of the rectangular barrier discussed in equa- 
tions (j9]l-( 111, the associated actions are elliptic, hyperbolic, or parabolic according to whether 



E is greater than, less than, or equal to Vo, respectively. 

To proceed further, let us note that by taking the conjugate of M„i with any matrix C e 
SU(1, 1); i.e., 

M«fe = CM„f,C-', (73) 
we get another matrix of the same type, forasmuch as Tr(M„i) = Tr(M„i). Conversely, if two 



transfer matrices have the same trace, one can always find a matrix C satisfying equation (73 1. 
The fixed points of M„f, are the image by C of the fixed points of Mq/,. In fact, if we write the 
matrix C as 

the matrix elements of Mq/, (marked by carets) and those of Maf, are related by 

(75) 



a|ci|2-ff*|c2|2-2/Im(;6cic*), 



P = y8c2-/?*c2-2/ciC2lm(Q'). 

In consequence, given any transfer matrix M„i one can always reduce it to a M„/, with one of the 
following canonical forms 



K(0) = 



exp(/<^/2) 
exp(-/0/2) 



A(^) = I '^^^^gg?,, '^^t/f/?? I' (76) 
^ ' sinh(^/2) cosh(^/2) ' 



N(y) = 



l-iv/2 v/2 
v/2 1+/V/2 



where < < 47r and ^, v e R. They have as fixed points the origin (elliptic), +/ and -/ (hyper- 
bolic) and +i (parabolic). All these SU(1, 1) matrices leave invariant - at each side 
of the potential, in agreement with equation (??). In addition, K(0) preserves the product X+X. 



A(^) preserves the quadratic form X\ and N(v) preserves the sum X+ + iX- ( ,Yonte et al. 
[2002| . 

The matrix K(0) represents the free propagation in a constant potential barrier with a dephas- 
ing of 0/2. Obviously, this reduces to a mere shift of the origin of phases. The second matrix 
A(^) represents a symmetric system with reflection and transmission phase shifts of = and 

= ±7r/2, and a transmission coefficient - sech(^/2). Finally, the third matrix, N(v), repre- 
sents a system having - cos(T|^)exp(iT|^) and - sin(T|^)exp(iT|;j), with tan(T|^) - v/2. There 
are many ways to implement these elementary actions depending on the physical system under 
consideration ( |Simon and Mukunda[ 11998"] ). 

The explicit construction of the family of matrices C is easy: it suffices to impose that C 
transforms the fixed points of M„/, into the ones of K((/i), A(f), or N(v). By way of example, let 
us examine the case when Mq^ is elliptic and its fixed point inside the unit disk is Zf- One should 
have 

O[CM„iC-',0] - ^[CUab,Zf] = <l>[C,z/] = 0. (77) 
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Rotation 
(Elliptic) 
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(Hyperbolic) 




(a) 

Parallel displacement 
(Parabolic) 




Figure 8: Plot of orbits in tlie unit disk for: (a) canonical transfer matrices in equation ^76) and (b) arbitrary transfer 
matrices. 



Solving this equation one gets directly 

1 



: exp(ij?) , 



C2 = -C1Z^, 



(78) 



and § is a real free parameter The same procedure apphes to matrices A(f) and N(y). 

The concept of orbit is especially appropriate for getting an intuitive picture of these actions. 
Given a point z, its orbit is the set of points z' obtained from z by the action of all the elements of 
the group. In figure [sja we have plotted typical orbits for each one of the canonical forms K(0), 
A(^), and N(v). They are 



z = C)[K((/i),z] =zexp(-/0), 

^rir.^ 1 ^-'tanh(^/2) 

z' = <I)[N(v),z] = 



1 +;ztanh(^/2) 
z + (l+iz)v/2 
l+(z-0v/2 ■ 



(79) 



For matrices K((p) the orbits are circumferences centered at the origin and there are no invariant 
hyperbolic lines. For A(^), they are arcs of circumference going from the point +/ to the point -/ 
through z and they are known as hypercicles. Every hypercicle is equidistant [in the sense of the 



distance ( 70 1] from the imaginary axis, which remains invariant (in the Euclidean plane the locus 
of a point at a constant distance from a fixed line is a pair of parallel lines). Finally, for N(v) the 
orbits are circumferences passing through the point +/ and joining the points z and -z* and they 
are denominated horocycles: they can be viewed as the locus of a point that is derived from the 



point +/ by a continuous parallel displacement (Coxeter 1968 1. 
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Figure 9: Unit liyperboloids defined in equation ^38^ , representing tlie space of states for SO(2,l). In eacli one of tliem 
we have plotted a typical orbit for the matrices A|^, A^, and A,^ (from left to right). In all the figures we have performed 
stereographic projection from the south pole 5 of the hyperboloid, to obtain the unit disk in the plane .v" = and the 
corresponding orbits, which are the actions of the SU(1, 1) matrices. 



For a general Mab the corresponding orbits can be obtained by transforming with the appro- 
priate matrix C the orbits deUneated before. In figure[8]b we have plotted examples of such orbits 
for elliptic, hyperbolic, and parabolic actions. We stress that once the fixed points of the transfer 
matrix are known, one can ensure that Za will lie in the orbit associated to Zb- 

An alternative way to understand these results is to look at the canonical matrices in the 



Lorentz group SO(2, 1). In fact, using (36 1 one finds that 



Af<(0) = 
A^(y) = 



1 
COS0 sin^ 
-sin0 cos<f. 

cosh^ -sinh^ 

1 
-sinh^ cosh^ 

( l+(v2/2) V -v2/2 

V I -V 



(80) 



l-(v2/2) j 



The action of these matrices in SO(2,l) is clear: A((^(ip) is a space rotation of angle (p in the 
1-2 plane, A^(^) is a boost in the direction of the axis 2 with velocity u - tanh^; and, finally, 
A|^(v) is a space rotation of angle T|;;| [such that tan(T|^) - v/2] followed by a boost of angle Tf^ 
and velocity v = tanh(v/2), both in the 1-2 plane. In figure |9] we have plotted examples of the 
orbits for each one of the subgroups in (80i. For A|^(0) the orbits are the intersection of the 
hyperboloid with planes = constant, for A^(^) with planes s' - constant, and for A^(v) with 
planes s" - = constant. Through stereographic projection we get the corresponding orbits for 
the matrices (|76ll in the unit disc. 



4.2. Transfer function in the half-plane 



The unitary matrix ( 16 1 plays an important role in the intertwining between the two basic 
vector bases used for the transfer-matrix description. In mathematical terms, 1A establishes a 
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one-to-one map between the groups SU(1, 1) and SL(2, R). To investigate the meaning of this 
map we observe that if the point vv e C is defined in terms of z by 



w = (D['W,z] = (81) 
\+iz 

then the interior of D is mapped onto the upper half-plane of the complex plane w, the boundary 
maps onto the real axis, while the exterior of D becomes the lower half-plane. This remarkable 
map is known as the Cayley transform. 
The metric now reads as 

ds^^- ^, (82) 

(Im wY 

and the geodesic lines are the open semicircles orthogonal to the real axis. The Mobius transfor- 
mations are 

bw + c 

Wh = <^{Mab, Wa\ = 7 , (83) 

hw + a 



with Mab obtained from M„/, through conjugation with U as in equation ( 17 1. They are also 
isometrics. 

The points w in the upper half-plane constitute the Poincare model of the hyperbolic plane 
H ( |Stahl|[T993[ ). In this way, one can transport all the geometrical properties of the unit disc D to 
the upper half-plane H. 

Since the matrix conjugation does not change the trace, the same geometrical classification 



in three basic actions still holds. In fact, by conjugating with t( the canonical forms (76 1, the 
corresponding ones for SL(2, R) are 



cos(0/2) sin(0/2) 
-sin(0/2) cos(0/2) 



e^l^ 





1 
V 1 



These matrices have as fixed points +i (elliptic), and oo (hyperbolic), and oo (parabolic), re- 
spectively. Clearly, "Ki^) is a rotation in phase space, also termed a fractional Fourier transform, 
while is sometimes called a squeezer or hyperbolic magnifier: it scales the positive amplitude 
+ up by the factor e^^^ and the negative one - down by the same factor. Finally, N(v) represents 



the action of a thin lens of power v (i.e., focal length 1/v) in geometrical optics (Wolf 20041. 



For the canonical forms ( 84 1, the orbits for a point w are 

cos(0/2)w - sin(<^/2) 



sin(0/2)vv + cos((p/2) 
w' = e'^w, (85) 



- w + v. 



In figure 10 a we have plotted these orbits. For matrices "^(0) they are circumferences cen- 
tered at the invariant point +i and passing through w and - 1 /w. For 'K'(^), they are fines going 

21 



Elliptic 



Hyperbolic Parabolic 





N 

N \ 












U • 


\ ', 








1 

1 














1 
I 





(b) 



Figure 1 0: Plot of typical orbits in the hyperbolic plane H: (a) canonical transfer matrices as in equation j84| and (b) 
arbitrary transfer matrices obtained by matrix conjugation. 

from to the oo through w and they corresponds to hypercicles. Finally, for matrices N(v) the 
orbits are lines parallel to the real axis passing through w and they are the horocycles (jCoxeter 



1969 1. It is in the plane H where the denomination of parallel displacements becomes manifest. 

As we did before, for a general matrix Mat the corresponding orbits can be obtained by 
transforming with the appropriate matrix C that transforms the fixed points of Mab into the ones 
of iCiip), ^(^), or Niv), respectively. In figure 10 b we have plotted typical examples of such 



orbits. 



4.3. Factoring the transfer matrix 

Many types of factorizations have been considered in the literature ( Arsenault and Macukow 
|1983[|Abe and Sheridan||1994{|Shamir and Cohen||1995| l, all of them decomposing the matrix as 
a unique product of other matrices of simpler interpretation. Particularly, given the relevant role 
played by the Iwasawa decomposition, both in fundamental studies and in applications to several 
fields, one is tempted to investigate also its role in the transfer-matrix formalism. 

Without embarking us in mathematical subtleties, the Iwasawa decomposition is established 
as follows ( jBarut and Rgczka} |1977[ IHelgason) |1978| l: any element of a (noncompact semi- 
simple) Lie group can be written as an ordered product of three elements, taken one each from 
a maximal compact subgroup K, a maximal Abelian subgroup A, and a maximal nilpotent sub- 
group A^. Furthermore, such a decomposition is global and essentially unique. 

For a matrix Mab £ SU(1, 1), the decomposition reads as 



M„, = K(,^')A(r)N(v'), 



(86) 



where the matrices appearing here are of the form of the canonical ones in equation ( |76| l, but the 
parameters (f>',^', and v' are given in terms of the elements a and /3 of the transfer matrix by 

<P'/2 = aTg(a + i/3), 
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Figure 1 1 : Representation in the unit disk of a transformation witli the parameters indicated in the text. The point Zb is 
transformed by the system into Za ■ The three orbits of the Iwasawa decomposition are indicated. 



r/2 = ln(l /\a + i/3\), (87) 
v'/2 = Re(a(3*)/\a + if3\^. 

The importance of the Iwasawa decomposition reflects at the geometrical level: no matter 
how complicated a system is, its action can always be viewed in terms of these three basic 
actions with a patent meaning. Let us show this with a practical example: we consider a system 



that transforms the point zi, - 0.4 exp(-/7r/3) into Za - -0.44 + 0.49 / [see Monzon et al. ( 2002 ) for 



a realistic implementation]. In figure 11 we have plotted these points Zb and Za in D. Obviously, 
from these data alone we cannot infer at all the path for this discrete transformation. 

The Iwasawa decomposition remedies this drawback: once we know the values of <p', 



and v' [that are easily computed from equation ( 87 1] we get the intermediate values of z' for 



the ordered application of the matrices K((f>'), Ai^'), and N(v'), which, in fact, ensures that the 



trajectory from zi, to z„ is defined through the corresponding orbits, as shown in figure 1 1 



4.4. Geometrical reflections as building blocks 

In the Euclidean plane any isometry is either a rotation or a translation. In any case, reflec- 
tions are the ultimate building blocks, since any isometry can be expressed as a composition of 
reflections. In this Euclidean plane, two distinct lines are either intersecting or parallel. Accord- 
ingly, the composition of two reflections in two intersecting lines forming an angle is a rotation 
of angle 20 around the intersection point, while the composition of two reflections in two parallel 
lines separated a distance c/ is a translation of value 2d. 

On the other hand, in hyperbolic geometry any two distinct lines are either intersecting (they 
cross in a point inside the unit disc), parallel (they meet at infinity; i.e., at a point on the boundary 
of the unit disc), or ultraparallel (they have no common points). A natural question arises: what 
is the composition of reflections in these three different kind of lines? To some extent, the 
answer could be expected: the composition is a rotation, a translation, or a parallel displacement, 
respectively. However, to gain further insights one needs to know how to deal with reflections in 
the unit disc. 

In the Euclidean plane, given any straight line and a point P which does not lie on the line, its 
reflected image P' is such that the line is equidistant from P and P' . In other words, a reflection 
is a special kind of isometry in which the invariant points consist of all the points on the line. 
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The concept of hyperbolic reflection is completely analogous: given the hyperbolic line f 
and a point P, to obtain its reflected image P' in £ we must drop a hyperbolic line £. from 
P perpendicular to I (such a hyperbolic line exists and it is unique) and extending an equal 



hyperbolic distance [according to ( 70 1] on the opposite side of XL from P. In the unit disc, this 
corresponds precisely to the notion of an inversion. 

We recall some facts about inversion ( |Coxeter| |1969[ ). Let C be a circle with center Q and 
radius R. An inversion on the circle C maps the point z into the point z' along the same radius in 
such a way that the product of distances from the center Q. satisfies 



k'-Q| k-Q|=7?2 



and hence 



Q + 



R^ + az* -Q*Q 
z*-Q.* ■ 



(88) 



(89) 



If the circle C is a hyperbolic line, it is orthogonal to the boundary of the unit disk and fulfills 
f2Q* -R^ + 1. In consequence 

.' = ^. (90) 

z - Ll* 



One can check (Pedoe 19701 that inversion maps circles and lines into circles and lines, and 
transforms angles into equal angles (although reversing the orientation). If a circle C passes 
through the points P and P' , inverse of P in the circle C, then C and C are perpendicular. More- 



over, the hyperbolic distance (70 1 is invariant under inversions. This confirms that inversions are 
indeed reflections and so they appear as the most basic isometrics of the unit disc. 

It will probe useful to introduce the conjugate bilinear transformation associated with a ma- 



trix Mab as [compare with equation ( 68 1] 



Za^^*{Mab,Zb\ = 



a*zl+(3* 
Pzl + a 



With this notation we can recast equation (90 1 as 



z' = 1>*[ln,z], 

where the matrix Iq e SU(1, 1) associated to the inversion is ( |Barriuso et al.[[2003bl l 

_/ -iSr iR i/R 
" " 1 -i/R i Q/R 



(91) 



(92) 



(93) 



The composition law for inversions can be stated as follows: if z' = <I'*[ln,z] and z" = <I'*[ln',z'] 
then 

z" ^^[\n'\n,z]. (94) 

To appreciate the physical meaning of the inversion, assume that incoming and outgoing am- 
plitudes are interchanged in the configuration shown in figure 1 . This is tantamount to reversing 
the time arrow. It is known that for a right-traveling mover X+, the conjugate amplitude [X+]* 
is a left phase-conjugate wave of the original one ( jZel'dovich et al. 1985| l. In other words, the 
time-reversal operation can be viewed in this context as the transformation 
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Elliptic 



Hyperbolic 



Parabolic 



Figure 12: Decomposition of tlie transfer-matrix action in terms of two reflections for tlie elliptic, hyperbolic and 
parabolic cases. 



for both a and h\ that is, it can be depicted by an inversion in the unit circle. The transformed 
points He outside the unit circle because time reversal transforms the upper sheet into the lower 
sheet of the hyperboloid. 

It is easy to convince oneself that the matrix relating these time-reversed amplitudes is pre- 
cisely M*^ and so the action can be put as 

i\ Za) ^ , , (96) 

which expresses a general property of the time-reversal invariance. 



In figure 12 we have drawn the three basic actions as a product of two reflections in two 
hyperbolic lines. For the elliptic case, the two hyperbolic lines Xi and X2 intersect at the fixed 
point Zf and form an angle 0, which is just one half of the rotation angle. The first inversion maps 
Zb into the intermediate point Zint, which is mapped into Za by the second inversion. Note that 
there are infinity pairs of lines satisfying these conditions, but chosen arbitrarily one of them, 
the other is uniquely determined. Once these lines are known, they delimit automatically the 
associated inversions. 

For the hyperbolic case, there are no invariant points in the unit disc, but the hyperbolic line t 
joining the fixed points Zf- and z/+ is the axis of the hyperbolic translation. We have also plotted 
the hypercicle passing through Zb and Za- The action can be now interpreted as the composition 
of two reflections in two ultraparallel hyperbolic lines Xi and X2 orthogonal to the translation 
axis. If Xi and X2 intersect the axis I at the points z\ and zi, they must fulfill 

d\i{zb,Za) = 2du(j.uZ7), (97) 

in complete analogy with what happens in the Euclidean plane. Once again, there are infinity 
pairs of lines fulfilling this condition. 

Finally, in the parabolic case, we have plotted the horocyle connecting Zb and Za and the fixed 
point. Now, we have the composition of two reflections in two parallel lines Xi and X2 that 
intersect at the fixed point z/, with the same constraints as before. 
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5. A closer look at the composition of transfer matrices 



5.1. Setting up the inverse system 



The property (19i has allowed us to characterize a compound system, as expressed more 



explicitly in (29 1. In its simplest form, it states that given two potentials V\ and V2, described by 
the transfer matrices M 1 and M2, with scattering amplitudes (r\ , t\ ) and (r2, f2), respectively (once 
more we drop the subscript ab to simplify the notation), the action of the compound system is 



M 



12 ■ 



M1M2 



and the reflection and transmission amplitudes associated to M12 are 
ri +r2exp(/2Ti) ^ tit2 



r\2 



1 +r|r2exp(i2Ti) 



hi 



1 +r*r2exp(/2Ti) 



(98) 



(99) 



with the same notation as in equation ( 44 1. With a bit of effort, one can derive valuable bounds 
on these coefficients ( [Visserl p[999 ; Boo nserm and Visser||2010a[ ). 

According to the general form (j7]i, the identity matrix has unit transmission and zero reflec- 
tion coefficients. In other words, it represents an antireflection system (without transmission 
phase shift). In consequence, two systems that are inverse, when composed give an antireflection 
system. 

Let us investigate the outstanding example of a single potential barrier of width L\ and height 
Vq. For the time being, we take E > Vq. Consequently, we look for another barrier of the same 
height such that, when put together with the original, gives the identity transfer matrix. 

A simple glance at the coefficients in (|9]), reveals two possibilities. The first is to couple other 
barrier with the same wave number k and width L2 such that 



sin[/(-(Li +L2)] = 0, 



(100) 



in such a way that the resulting barrier presents a transmission resonance. The second solution 
is to use a "complementary" barrier; that is, one having the same length and height, but opposite 
wave number -k, so it acts as canceling the effects of the first. 

Our geometrical picture leads to an appealing interpretation of these facts. In figure 13 we 
have schematized the action of the original barrier as a rotation of angle ip\ around the fixed 
points. The two previous solutions are a rotation of angle In-tpi and a rotation of angle -(fi, 
respectively, getting thus the identity in two different ways. 

The barrier can be seen as a quantum analog of a layer in classical optics. The condition 
£ > Vo ensures the classical regime in which light striking the layer is partly reflected and partly 
transmitted. On the contrary, E <Vo corresponds to an imaginary refractive index, producing 
total internal reflection (Bohm 1989[ l. So, the condition E >Vo prevents the appearance of total 
reflection. 

Negative values of the wave vector /c can be realized in terms of negative effective particle 
mass m ( |Kobayashil|2006[ [Dragoman and Dragoman||2007| l. The analog phenomenon in optics 
is a medium with both negative electrical permittivity e and magnetic permeability jj.. This is at 



the center of a lively and sometimes heated debate ( Cai and Shalaev 2009 1. This idea dates back 



to 1968, when Veselago (1968i theoretically predicted that these remarkable materials would 
exhibit a number of unusual effects derived from the fact that in them the vectors (k,E,B) of 
a plane wave form a left-handed (LH) rather than a right-handed (RH) set. For this reason, he 
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Figure 1 3 : Rotation in the unit disk associated witli a potential barrier described by the transfer matrix M i transforming 
the point zb = into Za = ri with the scattering amplitudes (ri = -0.3804 - i0.3339(, ti = 0.5689 - iO.6482). The orbits 
associated with the two solutions in the text can be either clockwise (of angle 2!T-<fi[ )oT counterclockwise (of angle ipi ) 
to give the identity. 
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Figure 14: Scheme of the energy flow for the system resulting by putting together two identical slab, one made of 
RH and other of LH material. Both constitute a pair of complementary media, each canceling the effect of the other. 
Consequently, no reflection occurs, as indicated by the dotted line. 



called them LH media. One of the most interesting properties of these LH materials is a negative 
refraction at the interface with a RH medium. Our solution for complementary barriers can be 



interpreted as putting together two RH and LH slabs (Mon zon et al.||2()()6||2008 |). The scheme 
of the energy flow in the resulting system appears in figure 14 



The discussion so far admits a straightforward generalization for any arbitrary potential. In- 
deed, let Miib be the transfer matrix of a potential that can be decomposed in an arbitrary number 
of barriers (some having positive and some negative values of the wave vector k), which can be 



constructed by a direct extension of ( 29 1. Now, we take the potential in the reverse order, which 
is represented by M/^, in (|8]l. Next, either we complete every barrier as in equation ( 100 1 or we 
switch every barrier with positive kj to an identical one with negative kj and viceversa. In both 
cases, this new system is represented by M^^. Since one can check that 



(101) 



when both systems are put together they give the identity. 
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This substitution kj —kj formalizes in a different framework the notion of "complemen- 



tary" media introduced by Pendry and Ramakrishna ( 2003| i: any medium can be optically can- 
celled by an equal thickness of material constructed to be an inverted mirror image of the 
medium, with e and fi reversed in sign. That is, complementary media cancel one another and 
become invisible (i.e., a perfect antireflector). 

5.2. The Wigner angle 

In special relativity there is an intriguing phenomenon that emerges in the composition of two 
noncollinear pure boosts: the combination of two such successive boosts cannot result in a pure 



boost, but renders an additional rotation, usually known as the Wigner rotation (Wyk 1984; Ben- 



Menaheml 1 19851 [Strandbergl [19861 |Aravmd| [T9971 [UngEI |200 1 [ [Malykinl |2006, .O 'Donnell an^ 



Visser 



1992^ 



iamilton 1996") is also used]. In other words, boosts are not a subgroup 



201 l|l [s ometimes the name of Thomas rotation (Jackson 1975 Ungar 1989 MuUer 



To fix the physical background, consider three reference frames K, K' and K" (see figure 15 1. 
Frames K-K' and K'-K" have parallel respective axes. Frame K" moves with uniform velocity 1)2 
with respect to K', which in turn moves with velocity ifi relative to K. The Lorentz transformation 
Ai2 that connects K with K" is given by the product L1L2, which can be decomposed as 



L1L2 = A12 = R(i/')L(i2)- 



(102) 



An equivalent decomposition in terms of a boost with the same modulus of iTn but with a diff'erent 
direction postmultiplied by the same rotation is also possible (Ungar 1989 Aravind 1997| l. 

In words, this means that an observer in K sees the axes of K" rotated relative to the ob- 
server's own axes by a Wigner rotation R(i/'). More explicitly, it is possible to show that the axis 



fi and angle of this rotation are ( ,Ben-Menahem |1985, ^Malykin. 2006, jRitus, |2008| ) 

?2 X 1^1 , . sin© 



|y2Xt?i| ' 

where is the angle between Ci and 1)2, and 



tan((A/2) : 



K + cos@' 



yi + 1 72 + 1 

yi - 1 72 - 1 



tanh2(^i/2)tanh2(^2/2)' 



(103) 



(104) 



yi and 72 being the corresponding factors for iJi and ij2, while and (2 are the rapidities. This 
means that tan(i/r/2) depends on the velocities as viU2, so the Wigner rotation is a second-order 
eff'ect and is absent in the non-relativistic limit. 

On the other hand, the resulting boost L(i2) has a velocity fulfilling 



ri2 =7172(1+^1 ■i'2) = 7172(1 +^1^2 COS©), 



(105) 



while the direction of ^12 has a complicated expression of little interest here. 

The set of boosts could be regarded as a hyperbolic space provided the velocity v is replaced 
by the hyperbolic parameter ( - tanh~'(i)), which constitutes the usual rapidity space and whose 
line element has a Lobachevskian metric, as known from long times ago (tLandau and Lifshitz 



|2000t |Rhodes and Semon|[2004| l. A triangle in this rapidity space obeys a non-Euclidean geom- 
etry and, in our context, this results in the fact that the parameters ( 105 1 of the compound boost 
L(i2) can be recast as 



cosh f 12 = cosh f 1 cosh (2 + sinh f 1 sinh ^2 cos , 
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Figure 15: (a) Reference frame K" moves with velocity !?2 relative to frame K' while frame K' moves with velocity 
Ci relative to frame K. (b) The axes of K" appears rotated relative to K by the Wigner angle 0. Time and one space 
dimension are suppressed for clarity. 



which is nothing but the hyperbolic law of cosines for the triangle induced by the boosts Li, L2, 
and L(i2). Therefore, given two sides and the included angle of the triangle (corresponding to the 
two non-coUinear boosts we wish to combine) one can determine the third side and its angle by 
a simple use of hyperbolic trigonometry. 

Moreover, a standard calculation shows that the expression 



103 I for the Wigner angle ij/ 



gives precisely the area of this triangle ( jChen and Ge 19981. We recall that for a hyperbolic 
(spherical) triangle the sum of the angles is less (greater) than n with the angular defect (excess) 
being the area. 

This suggests to look at the Wigner angle as a geometric phase. Roughly speaking, geo- 
metric phases are associated with the cyclic evolution of a system and the crucial concept to 
their understanding is anholonomy. Anholonomy (Shapere and Wilczek 1989 1 is a phenomenon 
in which non-integrability causes some variables to fail to return to their original values when 
others, which drive them, are altered round a cycle (the simplest example occurs in the paral- 
lel transport of vectors). This behaviour was anticipated by Pancharatnam when discussing the 
phase shift that appears in the coherent addition of two polarized beams on the Poincare sphere. 

It is worth mentioning that geometric phases associated with the group SU(1,1) have been 
previously identified ^Simon and Mukunda 1993a[ Monzon and Sanchez-Soto 1999 2001 1. The 
idea is to view the rapidity triangle as imbedded in the unit hyperboloid, which is a manifold of 
constant negative curvature (of value -1). The analogous triangle for rotations instead of boosts 
is traced on the unit sphere (of curvature -i-l) and the geometric phase appears as the area enclosed 
by the triangle on the sphere (Levi 1994 1. Thus, it is tempting to infer that the Wigner angle is 
just the area of the triangle on the hyperboloid, with the opposite sign to that of rotations: in fact, 
this is true as proved by|Aravind|(|1997|) and others pordan [|1988||Urbantkel |1990 1. 

After our discussion in the previous section, it is clear that the problem should arise in the 
context of transfer matrices. First of all, it is worth emphasizing that the two combining boosts 
and the resulting one are in the same plane, usually assumed for simplicity to be the 1-2 plane. 

In consequence, we restrict our attention to the composition of two Hermitian matrices H 1 
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Figure 16: Geodesic triangle in the unit disk for the potentials indicated in the text. By stereographic projection the 
associated triangle is also plotted in the unit hyperboloid. 



and H2 for, as explained before, they are the equivalent to pure boosts. In complete analogy with 
equation ( 



102 1 we have now 
H1H2 = M12 = UH(i2) 



exp(-/i/r/2) 








1 



exp(/i/'/2) / |fi2|\ ri2exp(-/i/r) 



1 



r* exp(/(/r) 



1 



where and f 12 are given by ( 99 1 and 



lA = 2argfi2 = arg(l H-rirj). 



(107) 



(108) 



The appearance of an extra unitary matrix is the signature of a Wigner rotation and, accordingly, 
the Wigner angle ijj viewed in SO(2,l) is just twice the phase of the transmission coefficient of 
the compound system. Obviously, when pi - p2, the Wigner rotation is absent, since then we are 
dealing with two parallel boosts, whose composition leads to the famous Einstein addition law 



of velocities (Vigoureux 1992 Vigoureux and Grossel 1993j). 



To show an explicit implementation of this phenomenon (|Monz6n and Sanchez-Soto 2001 



we take two potentials with Hermitian transfer matrices Hi and H2, and scattering coefficients 
in = 0.3736 - /0.2014,fi = 0.9055) and (r2 = 0.3413/,f2 = 0.9399). Equation ( fT07| fixes a third 
potential ^^^2) ^^'^^ when put together the compound system is an antireflection system with 



phase in transmission equal to the Wigner angle. This is shown in figure 16 where the triangle 
is also plotted in the unit hyperboloid. 

5.3. Hyperbolic turns 

According to Hamilton ( 1853| l, the turn associated with a rotation of axis n and angle )? is a 
directed arc of length )?/2 on the great circle orthogonal to n on the unit sphere. By means of 
these objects, the composition of rotations is described through a parallelogram-like law: if these 
turns are translated on the great circles until the head of the arc of the first rotation coincides with 
the tail of the arc of the second one, then the turn between the free tail and the head is associated 



with the resultant rotation ( jBiedenharn and Louck 1981 1. Hamilton turns are thus analogous for 
spherical geometry to sliding vectors in Euclidean geometry. It is unfortunate that this elegant 
idea of Hamilton is not as widely known as it rightly deserves. 
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Figure 17: Representation of the sliding turn Tf <-/2 in terms of two reflections in two lines Xi and X2 orthogonal to the 
axis of the translation I, which has two fixed points zf- and zf+- The transfomiation of a typical olf axis point zb is also 
shown. 



The purpose of this section is to show how the use of turns affords an intuitive and visual 
image of all problems involved in quantum scattering and reveals the emergence of hyperbolic 
geometry in the composition law of transfer matrices. 

Let us focus on the case of |Tr(M)| > 2. This is not a serious restriction, since any matrix of 



SU(1, 1) can be written (in many ways) as the product of two hyperbolic translations (Juarez and 
|Santander| [T9821 [Simon eTaL] |1989a| |Sanchez-Soto et al.| [2005 ; Sim on et al.||2006j l. The axis 
of the hyperbolic translation is the geodesic line joining the two fixed points. 

As explained in section 4.4, any pair of points z\ and zi on the axis of the translation I at 
a distance ^/2 can be chosen as intersections of Xi and X2 (orthogonal lines to €) with I. It is 
natural to associate to the translation an oriented segment of length ^/2, with 



(109) 



but otherwise free to slide on i (s ee figure |T7jl . This is analogous to Hamilton's turns, and will 
be called a hyperbolic turn T;_f/2 ( Simon et al. 1989b| l. 

Using this construction, an off-axis point such as Zh will be mapped by these two reflections 
(through an intermediate point Zint) to another point z„ along a curve equidistant to the axis. These 
other curves, unlike the axis of translation, are not hyperbolic lines. What matters is that once 
the turn is known, the transformation of every point in the unit disk is automatically established. 

Alternatively, we can formulate the concept of turn as follows. Let M be a hyperbolic trans- 
lation with Tr(M) positive [equivalently, Re(Q;) > 1]. Then, M is positive definite and one can 



ensure that its positive square root exists and reads as ( Barriuso et al. 2004 1 



VM : 



1 



VTr(M) + 2\ P* + l 



(110) 



This matrix has the same fixed points as M, but the translated distance is just half the induced by 
M; i.e., we set 

^(M) = 2^(VM). (Ill) 
This suggests that the matrix VM can be appropriately associated to the turn Te^^ji that represents 
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the translation induced by M. Therefore, we symbolically write 



Tf,f/2H^ VM. (112) 

Let ^1 and ^2 be the corresponding translated distances along intersecting axes ii and {2, 
respectively. Take now the associated turns Tfj_f|/2 and T;'2_^2/2 and slide them along ii and {2 
until they are "head to tail". Afterwards, the turn determined by the free tail and head is the turn 
associated to the resultant, which can be interpreted as a translation of parameter ^12 along the 
line (i2- 



This construction is illustrated in figure 18 where the pertinent parameters are (ri = 0.3103 - 



/0.8274,fi = 0.4383 + /0.1644) and (r2 = 0.6820 + /0.3079,f2 = 0.6601-/0.0659). The application 



of (|99| gives (rn = 0.5210 -/0.7331,fi2 = 0.3915-/0.1947). The noncommutative character is 
evident, and can also be inferred from the obvious fact that M12 M21. 

In Euclidean geometry, the resultant of this parallelogram law can be quantitatively deter- 
mined by a direct application of the cosine theorem. For any hyperbolic triangle with sides of 



lengths ^1 and ^2 that make an angle 0, the expression is precisely given in equation ( 106 1 



6. Periodic systems 

6.1. Finite periodic structures 

Periodic potentials are those whose shape is repeated indefinitely with period d; i.e., V(x) - 
V{x + d). The distinctive feature of these potentials is that the frequencies fall into continuous 
bands, separated by forbidden gaps. In the quantum context this was first noted by |Kronig and| 



Penney ( 1931 1 in the classic paper that laid the foundation for the modern theory of solids. This 
band structure also occurs, in principle, for mechanical, acoustical, and electromagnetic waves 
( [Griffiths and Steinke||200T] i. 

There is a recurring interest in the related instance where the potential V(x) consists of a 
finite number (say A^) of identical cells. We shall call that situation a finite periodic structure 
(the term locally periodic is also employed); such potentials are produced by any finite lattice 
and they are of great importance for a number of applications, such as superlattices, photonic 



crystals, multilayers, etc, where the finite size must unavoidably be taken into account (Felbacq 
|etaLl[T998i|Busch et al.||2007D . 




Figure 18: Composition of two hyperbolic turns Tf|^j/2 '^e-,,(i/2 by using a parallelogramlike law when the axes (i 
and £2 of the translations intersect. 
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From a theoretical standpoint finite periodic systems are more difficult to analyze because 



Bloch theorem ( jAshcroft and Mermin 1976 1, which so dramatically simplifies the periodic prob- 
lem, does not apply. It is amazing that the finite periodic case can be solved analytically for 
arbitrary A^. This was first discovered by Abeles (1948) and rediscovered in the quantum context 



by Kiang ( 1974) and|Cvetic and Pieman, (,1981 


) and later by several others ( 


Vezzetti and Cahay 


1986 Lee et aLj 1989 Kalotas and Lee |1991 


IGriffiths and Taussig 1992| 


Sprung et al. 


1993 


Wu et al. 1993 ; Rozman et al.V 1994b ; Liviotti 


1994; Chuprikov';i996i Erdos et al. 1997 


Barra 


and Gaspard,, 1999, .Sprung et al.,,1999,,Morozov et al.,,2002,,Pereyra and Castillo, ,2002 





Let us suppose that the arbitrary potential V(x) (the basic unit cell) is repUcated times at 



regular intervals, as schematized in figure 19 Our problem is to construct the transfer matrix for 
the whole array, given the transfer matrix M for the single cell. The amplitudes at the jth cell are 



Using this equation recursively we have that 

A-,0 



B-M 



(113) 



(114) 



so the whole problem reduces to the evaluation of M^. There are several elegant ways of cal- 
culating this power Perhaps, the most efficient is to use the Cayley-Hamilton theorem, which 



states that any square matrix satisfies its own characteristic equation (Gantmacher 2000| l. This 
means that 



M-'-luM + l =0. 



(115) 



where u - [Tr(M)] /2. Consequently, any higher power of M can be reduced to a linear combina- 
tion of M and the identity 1. By induction, we obtain the expression 



UN-l(u)M-UN-2(u)i. 



Here 



Un(0) ■ 



sm[(N+l)0] 
sin 6* 



(116) 



(117) 



with COS0 = M, are the Chebyshev polynomials of the second kind satisfying the recursion relation 
( Abramowitz and Stegun 1996| l 



U 



N+l 



^2uUn + U 



N>1, 



(118) 



and Uq(u) - I, Ui{u) = 2u. This provides a closed solution to the problem; in particular, for 
incidence from the left the reflectance = |rp of the array is 



[\/3\Un-i(u)V 



(119) 



which requires to know the transfer matrix for the unit cell. 

The Chebyshev polynomials play, for finite systems, a similar role to the one played by 
Bloch functions in the description of infinite periodic systems. In addition, they are very useful 
to perform numerical calculations. However, it is not easy to separate the different behaviors 
according to the value of the trace, as it happens for infinite periodic media. For this reason, we 
will follow an alternative route based on geometrical properties. 
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1 j N 

Figure 19: A finite periodic potential constructed from the basic cell V(x) 



6.2. Bandgaps in the unit disc 

First, we note that the reflectance associated to each one of the canonical matrices ( |76l ) is 



111 



tanh^(^/2), 



(120) 



= (y/2)2/[l+(v/2)2]. 



While is identically zero, and tend to unity when ^ and v, respectively, increase. 
Nevertheless, they have distinct growths: goes to unity exponentially, while "R^ goes as 



Let C be the matrix ( 74 1 that goes by conjugation from an arbitrary M to its canonical version. 
All the subgroups generated by K(0), A(^), or N(v) are one-parametric and therefore Abelian, so 
we have that 

U{^Ll)U{^i2)^U{^i^+^i2), (121) 

where // is the appropriate parameter (p, ^, or v. For an A^-cell system the overall transfer matrix 
is 

C"' C = C"' UiNn) C, (122) 

From this equation, one must expect three universal behaviors of the reflectance according the 
transfer matrix for the basic cell is elliptic, hyperbolic, or parabolic. We shall work in what 
follows the detailed structure of these three basic laws. 

Since the stop bands are given by the condition |Tr(M)| > 2, we first consider the case when 
M is hyperbolic. We can rewrite equation ( [73] l as 

M = C-'A(^)C, (123) 

and ^ is given by Tr(M) = 2 cosh^ > 2, because we are taking into account only positive values 



of Tr(M). One solution of equation ( 123 i is (Monzon et al. 2003 1 

ci = Fiji* + i sinh^) , C2 = -iF Im(a) , 



with F = 1/ V2sinh^[sinh^-Im(/3)]. 



Carrying out the matrix multiplications in (122 1 it is easy to compute the reflectance: 

|;S|2 + [sinh(^)/sinh(A?f)]2 ■ 
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(124) 



(125) 



Elliptic 



Hyperbolic Parabolic 



Figure 20: Successive iterates (N = 1, . . . , 5) for an elliptic, hyperbolic, and parabolic action starting from the origin as 
the initial point. The fixed points are marked in cyan. Only hyperbolic and parabolic actions tend to the unit circle. 



This is an exact expression for any value of A^. As grows, K^^^ approaches unity exponentially, 
as expected from a stop band. 

The band edges are determined by |Tr(M)| = 2; that is, when M is parabolic. A calculation 
very similar to the previous one shows that 

^('V) = ^ (126) 

^ I^P + (1/A^)2' ^ ^ 

with a typical behavior ' R!'^^ ~ 1 - 0{N~^) that is universal in the physics of reflection. The 



general results ( 125 i and ( 126 1 have been obtained in a different framework by ,Yeh, ( ,1988) and 
|Lekner| ( fT994l l. 

Finally, in the allowed bands we have |Tr(M)| < 2, M is elliptic, and 

^ Q'-2Qcos(2NE) 
l+<32-2<3cos(2A^H)' 



where 



e'==Re(a) + /Jl-[Re(a)]2. (128) 



Now the reflectance osciflates with between the values (Q^-2Q)/(Q - 1)^ and (tf + 2<3)/((3 + 

It seems quite pertinent to picture these behaviors in the unit disc. Note that if we have only 
an incident wave from the left (B- - 0; that is, Zh = 0) and simultaneously \za\ - 1 the system 
behaves as a perfect mirror Therefore, a mirror maps the origin into a point on the unit circle. 

Henceforth, we shall take zo = Zb = 0. As mentioned before, all the points zn obtained by 
iteration from zo lie in the orbit associated to the initial point zo by the single cell, which is 
determined by its fixed points. 



In figure 20 we have plotted the successive iterates worked out numerically for the three 
archetypical actions (Barriuso et al.|[2003a| l. In the elliptic case, the points zn revolve in the orbit 
centered at the fixed point and the system never reaches the unit circle. On the contrary, for the 
hyperbolic and parabolic actions the iterates converge to one of the fixed points on the unit circle. 



although with different laws, which correspond to the stop band and edges, respectively (Lekner 
[20001 ). 
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The Mh iterate can be easily computed for the canonical forms in equation ( 76 1 and then. 



conjugating as in (73 i. For a hyperbolic action one has 



zn = 



(129) 



i-eizfjzf-y 

where ^ - (a + /3zf-)/(a + /3zf+) is a complex number satisfying |^| < 1 and Zf± are the fixed 
points. Analogously, for the parabolic case we have 

,2 



Nf3z 



Zn - 



■f 



Npzf-l 



(130) 



Zf being the (double) fixed point. In both cases, zn converges to one of the fixed points on the 
unit circle, so \zn\ — > 1 when increases, a typical behavior of a stop band (or, in other terms, a 
perfect mirror). In the mathematical literature this limit point is referred to as the Denjoy- Wolff 
point of the map ( Kapeluszny et al. 1999| l. 



6.3. Bandgaps in the half-plane 

The previous formalism can be translated to the hyperbolic half -plane H by the unitary trans- 
formation ([16]), as explained in section 4.2. However, to give a physical feeling, we prefer to 
illustrate this point with the simple yet interesting case of optical beams ( Barriuso et al.||2005a| l. 

In paraxial-wave optics, axially symmetric (monochromatic scalar) beams are specified in 
the Hilbert space of complex-valued square-integrable wave-amplitude functions ^'{x) ( [Gloge 
and Marcuse 1969 1, with x labeling the axis. To deal with partially coherent beams we specify 
the field not by its amplitude, but by its cross-spectral density. The latter is defined in terms of 
the former as 



r(xi, ^2) = <^* (^1)^(^2)), 



(131) 



where the angular brackets denote ensemble averages. 

There is a wide family of beams, the Schell-model fields ("Wolf and Collett! '1978VFoley and| 
Zubairy|[T978l|Starikov and"WoIf| [19821 [Friberg and SudoT, 1982; Gori 1983; Gori and GreUaf 



1984 Friberg and Turunen 1988[(Ambrosini et al. 1994| l, for which the cross-spectral density 



( 1 3 1 1 factors in the form 



V{xi,X2)^ ^JI(xl)I(x2)fi(xl-X2). 



(132) 



Here / is the intensity distribution and jj. is the normalized degree of coherence, which is transla- 
tionally invariant. When these two fundamental quantities are Gaussians 



/(x) = 



I 



exp 



v2 ^ 



2crjj 



p(x) - exp 



v2 ^ 



2(7-2 y 



(133) 



the beam is said to be a Gaussian Schell model (GSM). Here, J is a constant independent of 
X that can be identified with the total irradiance and cr/ and cr^ are, respectively, the effective 
beam width and the transverse coherence length. Other well-known families of Gaussian fields 
are special cases of these GSM fields: when cr^ <K cr/ we have the Gaussian quasihomogeneous 
field, and when o-f^ ^ 00 we have the coherent Gaussian field. 

Anyhow, the crucial point for our purposes is that for GSM fields one can assign a complex 
parameter Q (Simon et al\ [T984[ [19^ [1988] [Simon and Mukundaj |1993b[ [Dragomanj [T996 
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Figure 21 : Plot of | Tr(A1)| in terms of the parameters gi and g2 of the optical resonator The plane | Tr(A^)| = 2 is also 
shown. The density plot of the three-dimensional figure appears at the top. 



where 




(135) 



and R is the wave front curvature radius. This parameter fully characterizes the beam and satisfies 
the Kogelnik abed law; namely, after propagation through a first-order optical system described 



by the matrix A1 as in ( 64 1, the parameter Q changes to Q' via 



e' = <D[M,e] = 



bQ + a' 



(136) 



On account of Im2 > by the definition ( 134i, one immediately checks that ImQ' > and we 



can thus view the action of the system as a bilinear transformation on the upper complex half- 
plane. When we use the metric ds - \dQ\/lmQ to measure distances, what we get again is the 
standard model of the hyperbolic plane H ( Stahl 1993 i. 



The whole real axis, which is the boundary of H, is also invariant under ( |136| l and represents 
wave fields with unlimited transverse irradiance (contrary to the notion of a beam). On the other 
hand, for the points in the imaginary axis we have an infinite wave front radius, which defines 
the corresponding beam waists. The origin renders a plane wave. 

The geometrical scenario presented before allows one to describe the evolution of a GSM 
beam by means of the associated orbits. Let us go back to the example of the optical cavity 
treated in section 3.3. The associated transfer matrix M fulfills 

TrM = 2(2^1^2-1). (137) 
Since the trace determines the fixed point and the orbits of the system, the g parameters establish 



uniquely the geometrical action of the cavity. To illustrate this point, in figure 21 we have plotted 
the value of | Tr(A1)| in terms of and g2- The plane | Tr(A1)| - 2, which delimits the boundary 
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Elliptic 



Hyperbolic 



Parabolic 



Figure 22: Successive iterates in the half plane H for typical elliptic, hyperbolic, and parabolic actions. For hyperbolic 
and parabolic actions, the iterates tend to the real axis. The fixed points are marked in cyan. 



between elliptic and hyperbolic action, is also shown. At the top of the figure, a density plot is 



presented, with the characteristic hyperbolic contours (Kogelnik and Li 1966 1. 



Assume now that the light bounces times through this system. The overall transfer matrix 
is M'^ and the transformed beam is represented by 



(138) 



where Qq is the initial point. 

Note that all the Qf^ lie in the orbit associated to 2o by the single round trip, which is 
determined by its fixed points. By varying the parameters g of the cavity we can choose to 
work in the elliptic, the hyperbolic, or the parabolic case (Baikal and Kim 2002| l. 

In figure 22 we have plotted the sequence of successive iterates for diff'erent kind of ray- 
transfer matrices, according to our previous classification. In the elliptic case, it is clear that the 
points Qf^ revolve in the orbit centered at the fixed point and the system never reaches the real 
axis. On the contrary, for the hyperbolic and parabolic cases the iterates converge to one of the 
fixed points on the real axis, although with different laws ( Barriuso et al.[[2005a| l. 

The iterates of hyperbolic and parabolic actions produce solutions fully unlimited, which are 
incompatible with our ideas of a beam. The only beam solutions are thus generated by elliptic 



actions and, according with equation ( 137 1, the stability criterion is 



0< 12^1^2-11 = I cos(0/2)|<L 



(139) 



where (p is the parameter in the canonical form "TC in equation ( 84 1. Such a condition is usually 
worked out in terms of algebraic arguments using ray-transfer matrices ( Siegman 1986[ l. 



6.4. Quasiperiodic sequences 

The most relevant property of periodic systems is the possibility of generating bandgaps 
for wave propagation. Spatial periodicity, however, is not necessarily a requirement for that: 
the presence of large bandgaps has been reported in aperiodic structures. Moreover, they share 
distinctive physical properties with both periodic (e.g., the formation of gaps) and disordered 
random media (e. g., the presence of highly localized states characterized by high-field enhance- 
ment and anomalous transport properties). 

In few words, aperiodic structures are made up of two (or more) incomplete periods stacked 
together to form an overall system that is neither fully periodic nor random, but somewhere in 
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between ("Macia 2009'). These structures can be roughly sorted into two classes, those that are 
quasiperiodic and those that are not. The former present a discrete Fourier spectrum character- 
ized by self-similar Bragg peaks (quite reminiscent of periodicity), whereas the latter usually 
exhibit complex diffraction spectra with singular scattering peaks, and multifractal scaling prop- 
erties. 

From a mathematical viewpoint quasiperiodic functions belong to the class of almost periodic 
functions (Bohr 19511. These functions can be uniformly approximated by Fourier series con- 



taining a countable infinity of pairwise incommensurate frequencies. When the set of frequencies 
required can be generated from a finite-dimensional basis, the resulting function is referred to as 
a quasiperiodic one. The simplest example of a quasiperiodic function is 



V{x) - cos X + cos{Ax) , 



(140) 



where A is an irrational number. 

Additionally, quasiperiodic structures lack translational invariance but possess a high degree 
of rotational symmetry, while the not quasiperiodic lack both translational and rotational sym- 
metry but display remarkable self-similarity (scale invariance symmetry) in their structural and 
spectral features ( |Macia| p006 ). 

The interest in these sequences was originally fueled by the theoretical predictions that they 



should manifest peculiar electron and phonon critical states ( [Ostlund and Pan dit 1984; Kohmoto 
et al.| |1987|l, associated with highly fragmented fractal energy spectra (Kohmoto et ai.| |1983 



T989t|Suto[|19"89l|Bemssard"eraLl[T989irOh and Lee|[TW3l|Chakrabarti et al.||1995TlLiu 



Luck 



1997 



Monsoriu et al. 2006 1. On the other hand, the practical fabrication of Fibonacci (Merlin 
1985 ) and Thue -Morse (Me rlin et al.||198T ) superlattices triggered a number of experi- 
mental achievements that have provided new insights into the capabilities of quasiperiodic struc- 



et al. 



tures (Zarate and Velasco 2001 Velasco and Garcia-Moliner 2003 [ Albuquerque and Cottam 



2003 1. In particular, possible optical applications have deserved major attention and some in- 



triguing properties have been demonstrated (Tamura and Nori!ll989llAvishai and Berend 


1990 


Kolaf et al. 


19911 


Hattori et al. 1994 Vasconcelos and Albuquerque 1999; Hollingworth et al. 


2001 Lusk et al.| 


2001 Barriuso et al. 


2005bt |Mugassabi and Vourdas 


2009|l. Underlying all 



these theoretical and experimental efforts a crucial fundamental question remains concerning 
whether quasiperiodic devices would achieve better performance than usual periodic ones for 
some specific applications ( Macia 2001 Barriuso et al.||20()3a i. 

A simple understanding of a well ordered but aperiodic arrangement of numbers can be 
grasped by thinking of the Fibonacci numbers sequence F„ = {1, 1,2,3,5,8, 13,21, .. .). The terms 
in this sequence are generated from the recursive equation F„+i = F„ +F„-i, starting with Fq-I 
and Fi = I. Hence, each number in the sequence is just the sum of the preceding two. The 
sequence is then perfectly ordered, but the rule used to generate it has nothing to do with period- 
icity. The symbolical analog of the Fibonacci sequence, constructed using two types of building 
blocks, say A and B, can be obtained from the substitution rule A i-> AB and Z? i-> A, whose suc- 
cessive application generates the sequence of letters A,AB,ABA,ABAAB,ABAABABA, . . . and so 
on. In this way, we get a perfectly ordered word which is not periodic at all. 

The standard method of constructing aperiodic structures is thus through a substitution rule 
operating on a finite alphabet {A,B, . . .}. which consists of certain number of letters. In actual 
realizations each letter will correspond to a different type of building block in the structure. In 
particular, the substitutional sequences that act upon a two-letter alphabet {A,B} are especially 
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Table 1 : List of the substitution rules determining the sequences usually considered in the study of self-similar aperiodic 
systems. 



Sequence Alphabet Substitution rule 



Fibonacci 


{A,B} 


o-i 


^AB, 0-2 = 


A 


Thue-Morse 


{A,B] 


o-i 


= AB, 0-2 = 


BA 


Period doubling 


{A,B} 


o-i 


^AB, 0-2 = 


AA 


Rudin-Shapiro 


{A,B,C,D] 


o"i 


= AC, 0-2 = 


DC, 0-3 = AB, 0-4 = DB 


Circular 


{A,B,C} 




= CAC, 0-2 


= ACCAC, 0-3 = ABCAC 



important: in this case the algorithm reduces to 



Ai-^a-i{A,B), 



Bh^a-2{A,B), 



(141) 



where a-\ and 0-2 can be any string of A and B. In table[T]we list some representatives among the 
plethora of aperiodic structures grown during the last two decades. 

To each rule we can associate a substitution matrix T, whose columns give the number the 
letters A and B which occur in the substitutions cr\ and 0-2 



nA[cry{A,B)] nA[cr2(A,B)] ^ 
nB{(ri{A,B)] nB{cT2iA,B)] 



(142) 



This matrix does not depend on the precise form of the substitutions (the order of the letters), 
only on the number of letters A or B. 

The eigenvalues of the substitution matrix T contain a lot of information. In fact, according 
to a theorem by Bo mbieri and TayIorj ( |1986) l, if the spectrum of T contains a Pisot number, the 
structure is quasiperiodic; otherwise it is not. A Pisot number is a positive algebraic number (i.e., 
a number that is a solution of an algebraic equation) greater than one, all of whose conjugate 
elements (the other solutions of the algebraic equation) have modulus less than unity (Bertin| 

1992| l. For example, let us consider the sequence crj - AAAB, 



et al. 1992 Godreche and Luck 



1X2 = BBA, whose characteristic matrix is: 



3 1 
1 2 



(143) 



The eigenvalues are /li = (5 + V5/2) and A2-{5- V5/2). Since A\ > /I2 > 1, the sequence does 
posses the Pisot property. In Table |2] we give a sketch of the properties of the substitutional 
sequences considered thus far 

The sequences can be also characterized by the nature of their Fourier spectrum ( |Severin| 
1989| l. The Fourier spectrum corresponding to a perfect infinite periodic system 



and Riklund 



contains delta functions centered in wave numbers associated to the reciprocal lattice (this is the 
origin of the Bragg peaks). On the contrary, a disordered structure has a very flat spectrum. Ape- 
riodic heterostructures following a deterministic sequence display characteristic spectral proper- 
ties absent in either of these extreme cases. 

For a specific sequence of length A^, the discrete Fourier transform is 



N-\ 



1 \— ' / —2mk 
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Sequence Substitution matrix Eig 


envalues 


Fibonacci 


' 1 








^1 


= (1+ V5)/2, ^2 = (1- V5)/2 


Thue-Morse 


' 1 


1 ^ 






^1 


= 2, i2 = 


Period doubling 


,2 


) 






^1 


= 2, i2 = -l 




f 1 





1 


^ 






Rudin-Shapiro 



1 




1 


1 




1 





= 0, ^2 = 2, /I3 = V2, /I4 = - V2 




^ 


1 





1 J 








^ 1 





2 ^ 








Circular 


2 





3 




^1 


= -1,^2 = 2+ V(5), /l3=2- V5 




^ 2 


1 


2 J 









where w(y') is a numerical sequence obtained by assigning to each letter of the alphabet a fixed 
amplitude. This assignment is otherwise arbitrary and does not change any conclusion. In con- 
sequence, one could, e.g., use A i-> -1 and B i-> 1. The structure factor (or power spectrum) is 
( [Cheng and Slyit] [1990) 1 



(145) 



From a rigorous viewpoint, the only well-established concept attached to the Fourier spec- 
trum is its spectral measure. If we define 



dvfiik) - S N(k)dk, 
we will be concerned with the nature of the limit 



dv(k) - Hm dvN(k), 



(146) 



(147) 



which corresponds to an infinite structure and a continuous variable k. Just as any positive 



measure, ( 147 1 has a unique decomposition (Reed and Simon 1980j ) 

dv(k) - dvpp(k) + dy^cik) + dv^dk) 



(148) 



into its pure point, absolutely continuous and singular continuous parts. 

The pure point part refers to the presence of Bragg peaks; the absolute continuous part is a 
difFerentiable function (difiiise scattering), while the singular continuous part it is neither con- 
tinuous nor does it have Bragg peaks. It shows broad peaks, which are never isolated and, with 
increasing resolution, split again into further broad. 

The Fibonacci sequence has a pure point spectrum; the Thue-Morse sequence has a singular 
continuous Fourier spectrum, while the Rudin-Shapiro sequence shows an absolute continuous 
one. For a very detailed and up-to-date discussion of these issues, the reader is referred to 
comprehensive book by Macia ( 2009| l. 
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6.5. Hyperbolic tilings 

The rich properties of these aperiodic structures suggest the utility of studying of systems 
based on more general sequences ( jSpinadef) |1999| l. Periodicity is intimately connected with 
tessellations, i.e., tilings by identical repUcas of a unit cell (or fundamental domain) that fill the 
plane with no overlaps and no gaps. Of special interest is the case when the primitive cell is a 



regular polygon with a finite area ( Zieschang et al. 1980 1. In the Euclidean plane, the associated 
regular tessellation is generically noted {p, q], where p is the number of polygon edges and q is 
the number of polygons that meet at a vertex. Geometrical constraints limit the possible regular 
tiUngs {p,q] to those verifying 

(p-2)(q-2)^4. (149) 

This includes the classical tilings {4,4) (tiling by squares) and {6,3) (tiling by hexagons), plus a 
third one, the tiling {3,6) by triangles (which is dual to the {6,3)). 

On the contrary, in the hyperbolic disk regular tilings exist provided (p - 2){q - 2) > 4, which 
now leads to an infinite number of possibilities (Magnus 1974| l. The fundamental polygons are 
connected to the discrete subgroups of isometries (or congruent mappings); they are Fuchsian 
groups (|Ford| 1 1972)1 and play for the hyperbolic geometry a role similar to that of crystallographic 



groups for the Euclidean geometry (Beardon 1983 1 



A tessellation of the hyperbolic plane by regular polygons has a symmetry group that is 
generated by reflections in geodesies, which are inversions across circles in the unit disc. These 
geodesies correspond to edges or axes of symmetry of the polygons. Therefore, to construct a 
tessellation of the unit disk one just has to built one tile and to duplicate it by using reflections in 
the edges. 

To go straight to the point let us consider the following parabolic transformations 



1-/ 
1 



1 

1+/ 



1+/ 
1 



1 

1-/ 



(150) 



with fixed points +/ and -/, respectively. A possible implementation of these matrices (and their 
inverses) in terms of two commonly employed materials can be found in Barriuso et al. (2009 1. 
In figure 23 we have plotted the tessellation obtained by transforming the fundamental square 
with the Fuchsian group generated by the powers of {A,B} (and the inverses). 

To give an explicit construction rule for the admissible words, we proceed as follows. First, 
we arbitrarily choose a side of the fundamental square and assign to it the value 0. Then the other 
three sides are numbered clockwise as 1, 2 and 3. It is easy to convince oneself that this assign- 
ment fixes once and for all the numbering for the sides of all the other squares in the tessellation. 
However, these squares can be distinguished by their orientation (as seen from the corresponding 
center): the clockwise oriented ones are filled in red, while the counterclockwise ones are filled 
in yellow. In short, we have determined a fundamental coloring of the tessellation ( |Griinbaurii| 
and Shepard|[T987l l. 



To derive a center Zn+i from a previous one z„, one looks first at the corresponding color 
jump. Next, the matrix that takes z„ into z„+i depends on the numbering of the side (0, 1, 2, or 3) 
one must cross, and appears in the column labeled T in table [3] The next generation is obtained 
in much the same way, except for the fact that A„ and B„ must be replaced by A„+i and B„+i , 
respectively, as indicated in the table. In creating recursively any word, the origin is denoted as 
Zo and the matrices Aq and Bq coincide with A and B. One can then construct any word step by 
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Figure 23: Tessellation of the unit disk with the matrices jl50) . The marked points are the centers of the squares. All of 
them are the transformed of the origin by a matrix that have as reflection coefficient the complex number that links the 
origin with the center of the square. 



Table 3: Explicit rules to obtain the center z.„+[ from z„ in the tiling by hyperbolic squares. We have indicated the 
corresponding transformations, which depend on the color jump and the sides crossed by going from z„ to z„+i . 



Side 


T 


red yellow 


T 


yellow — i 
A„+i 


red 

Bn+l 





A„ 


B„ 




A-' 


An 


An B„An 


1 


B„ 


B„A„ B~' 


B„ 


B-1 

n 


b,7'a„b„ 


B„ 


2 


B-J 


B„ ^A„B„ 


B„ 


B„ 


B„A„B„ ' 


B„ 


3 




A„ 


A„ B„A„ 


An 


An 


AnBnA„ 



Step. For example, the word that transforms zo into ze in the zig-zag path sketched in figure [23] 
results; 
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-^zi 


A, 








Zl 


-^Z2 


AB" 


'A- 






Z2 


-^Z3 


AB- 


'A- 






Z3 


-^Z4 


AB- 


'B- 


ABBA 


-1 


Z4 


-^Z5 


AB- 


'B- 


ABBA 


-1 ' 


Z5 


->Z6 


AB- 


'B- 


AAB-i 


A-'A-^BBA 



The characteristic matrices of these substitution rules are 



I odd ■ 



( 1 





1 


1 ^ 


2 


1 


3 


1 


1 


1 


1 





. 3 


1 


2 


1 J 



1 1 








^ 





1 














1 





. 








1 J 



(152) 



with eigenvalues 2 + VS, 2 - VS, 1 and - 1 for Todd and 1 for Tgven (that produces a trivial effect 
and will no be taken into account). Since one of the eigenvalues is greater than 1 and the others 
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Figure 24: Normalized power spectrum for a word in the zig-zag patli sliown in figure[23]but witli 1600 letters. 



have an absolute value less than or equal to unity (with at least one of modulus 1), the substitution 



possesses the Salem property (Salem 1963 i, property weaker than the Pisot property. 



In figure 24 we have plotted the structure factor of this sequence for a word of 1600 letters 
that starts at the origin. Previously, we assigned to each letter in the alphabet the quartic roots of 
the identity, namely A i-> /, B i-> +1, A~' i-> -/, B~' i-> -1. The gross features of the spectrum 
are seen to be humps separated by almost empty regions. Inside these humps, there is a blurred 
structure built up of packed delta-spikes. The dominant peaks tend to be isolated and larger 

To describe the continuous part of the measure in an analytical way seems to be a hard task. 
Still we can introduce the integrals ( jAubry et al. 1987 Godreche and Luck 1990| l 



Jo 



ik)dk. 



(153) 



where kmax is some arbitrary spectral cutoff. The nature of the limit measure is partly coded in 
the behavior of these quantities for large A^. In fact, for an absolutely cont inuous measure they 
tend to be constant, while they scale as 1^ ~ N~^^^ for a pure point measure ( Godreche and Luck 



19921. A somehow intermediate behavior can be expected for the singular continuous case. For 



our system, we have evaluated numerically ( 153 1, finding the law 



7 = 0.36 ±0.02. 



(154) 



This rules out the existence of an absolutely continuous component and suggests that only a 
singular spectrum is present. This exponent can also be related to the theory of multifractals: if 
the measure dv(k:) has a generahzed dimension function D^, then y should be linearly related to 
Dk Pentschel and Procaccial|1983| l. 



7. Concluding remarks 

This review is concerned with the transfer matrix, a powerful tool that relies only on the 
linearity of a system with two input and two output channels. Therefore, it is not surprising the 
variety of domains in which this object has been successfully employed. 
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Instead of embarking in a detailed description of a particular model, our thread has been to 
put geometry to work. This provides a useful and, at the same time, simple language in which 
numerous physical ideas may be clearly formulated and effectively treated. 

In a first step, we have transplanted the transfer matrix into space-time phenomena. This 
gateway works in both directions: here, it has allowed us to establish a relativistic presentation 
of the transfer matrix, but specific models can be also used as an instrument for visualizing 
special relativity. This is more than an academic curiosity: in fact, some intricate relativistic 
eflTects, such as, e.g., the Wigner angle, can be measured (and not merely inferred) by optical 
setups. 

By resorting to elementary notions of hyperbolic geometry, we have interpreted in a natural 
way the action of the transfer matrix as a mapping on the unit disk and on the upper half-plane. 
The trace turns out to classify and characterize the basic geometrical actions, which has physical 
relevance. In fact, in this arena, nontrivial phenomena can be understood in terms of analogues 
of vectors in Euclidean geometry. 

We have applied this perspective to periodic systems, explaining the existence of bandgaps 
in appealing terms. In addition, we have presented schemes to generate quasiperiodic sequences 
based on tessellations of the unit disc. 

Nothing of the material presented here is applicable per se, but everything can be relevant for 
researchers in other fields. This is the beauty of the approach. To our mind, the manuscript is 
better than ever. May each reader benefit and enjoy! 
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